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— '?1^-^o0LY-ffisT0et,"l2ations  at  the  KC‘SCF/CI  Levet 

I.  This  past  year  we\enlarged  for  the  CRAY  XMP  (COS  and  CTSS)  the  MRD-CI  program 
(to  include  more  activ?  orbitals,  more  reference  conf igurations ,  more  roots  and 
llarger  basis  sets. /An  MRD-CI  population  analysis  program  was  written.  A  program 
Iwas  written  tojodate  centroids  of  localized  molecular  orbitals.  We  implemented 
and  initiate«Kprototype  calculations  on  a  major  program  system  GAMESS  which  has 
the  capability  of  optimizing  molecular  geometries  by  a  derivative  method  at  the 
SCF  MC-SCF/CI  and  CASSCF  levels. 

] developed  and  implemented  a  new  approach  based  on  localized  occupied  and 
virtual  orbitals  in  the  interaction  region  with  the  remainder  of  the  localized 
molecular  orbitals  being  folded  into  an  effective  Cl- Hamiltonian.  This  method  is 
completely  general  and  applicable  to  reactions  and  molecular  decompositions  of 
energetic  compounds  for  ground  and  electronically  excited  states. 

We  initiated  a  major  research  effort  on  ab-initio  MRD-CI  calculations  on 
opening  the  ring  of  oxetane  or  protonated  oxetane  and  on  the  opening  of  the 
protonated  oxetane  ring  in  the  course  of  interaction  with  oxetane. 

We  have  carried  out  MRD-CI  calculations  for  a  great  many  geometries  for 
various  paths  of  attack  for  the  system  oxetane  plus  protonated  oxetane.  We  have 
varied  the  angle  of  opening  of  the  protonated  oxetane  ring  (ring  'A'),  varied  the 
internuclear  distance  between  the  rings  and  the  angle  between  the  plane  of  the 

oxetane  (ring  ‘ B ’ )  and  the  plane  of  the  protonated  oxetane  ring./'" _ Q 
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III.  The  0  NMR  spectra  of  several  oxetanes  were  examined  experimentally  by 
Walter  S.  Koski  to  see  if  they  could  be  used  to  calibrate  our  computed  electros¬ 
tatic  molecular  potential  contour  (EMPC)  maps  for  the  basicity  of  oxetane  and  its 
derivatives. 
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IV.  We  recently  had  the  opportunity  to  validate  the  reliability  both  of  the  SCF 
intermolecular  interaction  calculations  usin^  our  ab-initio  MODPOT  basis  set  and 
of  the  semi-theoretical  expressions  we  had  been  using  for  dispersion  energy  by 
comparing  our  results  on  the  nitromethane  dimer  with  those  of 

recent  larger  basis  set  [4,3, 1/2,1]  all-electron  calculations  on  the  nitromethane 

dimer,  [SCF,  SCF/SDQ-MBPT(4) ,  SCF/di spersion  energy  E^^]. 

Our  results  further  suggest  that  it  appears  feasible  to  correct  for  the  small 
differences  in  the  absolute  intermolecular  interaction  energies  between  the  ab- 
initio  MODPOT  SCF  results  and  those  of  all-electron  larger  basis  set  SCF 

calculations  by  the  difference  in  the  calculations  with  the  different 

basis  sets.  In  addition,  the  semi-theoretical  method  we  have  been  using  to  es¬ 
timate  dispersion  energies  proved  to  give  values  very  close  to  Egjj^  (variation- 
perturbation)  that  had  been  calculated  using  a  larger  all-electron  basis  set. 


V.  We  explored  a  number  of  featuresof  the  GAMESS  program  which  has  the  capability 
of  optimizing  geometries  by  a  derivative  method  at  the  MC-SCF  and  CASSCF  levels  as 
well  as  at  the  SCF  level.  We  are  exploring  and  using  more  of  the  new  options  in 
the  GAMESS  program  for  our  studies  on  the  energetic  oxetanes. 

VI.  Commensurate  with  the  ONR  priorities  expressed  to  us  by  our  ONR  Contract 
Monitor  to  devote  our  major  attention  to  the  problem  of  the  cationic  polymeriza¬ 
tion  of  oxetanes,  we  only  devoted  very  minimal  effort  to  further  development  and 
testing  of  the  POLY-CRYST  program  (for  ab-initio  calculations  of  crystals  and 
pol ymers ) . 

For  efficiency  we  initiated  changes  in  the  way  the  integrals  are  picked  up. 

We  also  ran  a  few  test  calculations. 
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CONCISE  SUMMARY 


I.  Conversion  of  Additional  Programs  to  CRAY  Vector  Supercomputers  and 
Program  Enhancements 

Our  exploratory  ab-initio  MRD-CI  (multireference  double  excitation- 
configuration  interaction)  calculations  on  the  reaction  of  oxetane  with 
protonated  oxetane  for  the  propagation  step  in  cationic  polymerization 
indicated  that  even  for  oxetane  and  certainly  for  the  substituted  oxetanes 
we  were  going  to  have  to  consider  more  active  orbitals,  more  reference 
configurations,  more  roots  and  larger  basis  sets  than  the  MRD-CI  program  was 
capable  of  handling  at  that  time. 

Thus,  this  past  year  we  enlarged  these  aspects  of  the  MRD-CI  program 
for  the  CRAY  XMP  Computers  (for  both  the  COS  and  CTSS  operating  systems). 

A  population  analysis  program  to  calculate  charges  on  the  atoms  and 
total  overlap  populations  between  atoms  from  the  density  matrix  following 
the  final  Cl  calculations  was  also  written. 

To  help  in  choosing  molecular  orbitals  most  optimally  localized  in  the 
reaction  (bond  making  and  bond  breaking)  region,  we  wrote  a  program  to 
locate  the  centroids  of  each  localized  molecular  orbital. 

This  past  year  we  implemented  and  initiated  prototype  calculations  on  a 
major  program  system  GAMESS  (Gaussian  Atomic  and  Molecular  Energy  Systems) 
which  has  the  capability  of  optimizing  geometries  by  a  derivative  method  not 
only  at  the  SCF  level  but  also  at  the  MC-SCF/CI  (multi-configuration 
SCF/conf iguration  interaction)  or  CASSCF  (Complete  Active  Space  MC-SCF) 
levels.  This  can  prove  important  in  any  calculations  on  energetic  compounds 
involving  nitro,  nitrato  and  nitramino  groups. 

As  a  result  of  our  expertise  on  CRAY  computers,  we  acted  as  "friendly 
users"  on  the  NSF  CRAY  XMP-48  supercomputer  at  the  San  Diego  Supercomputer 
Center  (SDSC)  located  at  GTA  San  Diego,  California  for  December  1985  prior 
to  the  official  opening  of  that  center. 

As  a  result  of  our  interaction  with  that  SDSC  computer  center,  by 
invitation  we  have  been  serving  on  their  review  committee  for  granting  NSF 
computer  time  on  their  CRAY  XMP-48. 


II.  MRD-CI  Calculations  for  the  Propagation  Step  in  Cationic  Polymerization 

of  Energetic  Oxetanes 

Our  major  research  this  past  year  has  been  on  theoretical  quantum 
chemical  ab-initio  MRD-CI  (mul tireference  double  excitation  -  configuration 
interaction)  calculations  on  the  mechanism  of  the  propagation  step  of 
cationic  polymerization  of  oxetane  by  protonated  oxetanes.  The  oxetane  - 
protonated  oxetane  systems  (especially  the  substituted  ones)  are  larger  than 
any  molecular  reacting  systems  for  which  any  MRD-CI  calculations  have  been 
carried  out.  We  developed  and  implemented  a  new  approach  based  on  localized 
occupied  and  virtual  orbitals  in  the  interaction  region  with  the  remainder 
of  the  localized  molecular  orbitals  being  folded  into  an  effective  Cl 
Hamiltonian.  This  method  is  completely  general  and  applicable  to  reactions 
and  molecular  decompositions  of  energetic  compounds  for  ground  and 
electronically  excited  states. 


This  past  year  we  initiated  a  major  research  effort  on  ab-initio  MRD-CI 
calculations  on  opening  the  ring  of  oxetane  or  protonated  oxetane  and  on  the 
opening  of  the  protonated  oxetane  ring  in  the  course  of  interaction  with 
oxetane.  We  considered  the  following  localized  orbitals  explicitly  in  the 
MRD-CI  calculations:  in  protonated  oxetane  'A'  ring  C^-0^,  ^ A~^1A’  the 


lone  pair  on  0 
in  oxetane 


1A’  ”1A 
'  B '  ring,  0 


01fl-H+,  and  the  bonds  connecting  hydrogens  to  C^A  and  C 


'4A 


0 1 B  C4A. 


1B~  zB ’ 


0 1 8_C4B ' 


the  lone  pair  on 


01B  and 


inter-ring 


We  have  carried  out  MRD-CI  calculations  for  a  great  many  geometries  for 
various  paths  of  attack  for  the  system  oxetane  plus  protonated  oxetane.  We 
have  varied  the  angle  of  opening  of  the  protonated  oxetane  ring  (ring  'A'), 
varied  the  internuclear  distance  between  the  rings  and  the  angle  between  the 
plane  of  the  oxetane  ring  (ring  * B ' )  and  the  plane  of  the  protonated  oxetane 
ring. 


The  most  favorable  pathway  we  have  found  is  for  to  attack  linearly 
along  the  C^-0^  bond  direction  in  an  type  of  attack  with  at  least  a 
partial  inversion  of  the  H  atoms  on  C^  with  concomitant  opening  of  the 


C4A~CIA  bond- 


III, 
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0  NMR  Spectra  of  Oxetanes 


The  0  NMR  spectra  of  several  oxetanes  were  examined  experimental ly  by 
Walter  S.  Koski  to  see  if  they  could  be  used  to  calibrate  our  computed 
electrostatic  molecular  potential  contour  (EMPC)  maps  for  the  basicity  of 

oxetane  and  its  derivatives.  The  results  for  oxetane  and  BAMO  show  that  ^0 
NMR  chemical  shifts  cannot  be  used  directly  to  make  such  calibrations.  The 
reasons  for  this  is  that  two  opposing  effects  contribute  to  the  resultant 
chemical  shift.  One  is  a  diamagnetic  effect,  related  to  the  electron 
density  on  the  0  nucleus,  and  a  second  is  a  paramagnetic  effect,  due  to  the 


w 


a 


presence  of  low  lying  excited  states.  It  is  difficult  to  unscramble  these 
two  contributions  and  get  a  direct  measure  of  the  electron  population  on  the 
oxygen  atom.  The  sign  of  the  paramagnetic  term  is  opposite  to  that  of  the 
diamagnetic  effect.  The  effects  of  these  two  terms  are  difficult  to 
disentangle  quantitatively. 

17  13 

However,  a  detailed  0  and  C  NMR  study  of  substituted  oxetanes 
coupled  with  other  chemicophysical/physicochemical  measurements  would  give 
us  a  great  deal  of  insight  into  the  electronic  properties  of  the  series  of 
energetic  oxetane  derivatives  of  interest  to  this  ONR  program. 


IV.  Ab-initio  Atom-Class  -  AtGm-Class  Potential  Functions  and  Independent 
Confirmations  of  the  Validity  of  Our  Ab-Initio  MODPOT/VRDDO  Energy 
Partitioned  Method  For  Intermolecular  Interactions. 

We  recently  had  the  opportunity  to  validate  the  reliability  both  of  the 
SCF  intermolecular  interaction  calculations  using  our  ab-initio  MODPOT  basis 
set  and  of  the  semi-theoretical  expressions  we  had  been  using  for  dispersion 
energy  by  comparing  our  results  on  the  nitromethane  dimer  with  those  of 
recent  larger  basis  set  [4, 3, 1/2,1]  all-electron  calculations  on  the 

nitromethane  dimer,  [SCF,  SCF/S0Q-MBPT(4) ,  SCF/di  spersion  energy  E^?^  ]  . 

Our  results  further  suggest  that  it  appears  feasible  to  correct  for  the 
small  differences  in  the  absolute  intermolecular  interaction  energies 
between  the  ab-initio  MOOPUT  SCF  results  and  those  of  all-electron  larger 

basis  set  SCF  calculations  by  the  difference  in  the  E^MTP  calculations 

with  the  different  basis  sets. 

In  addition,  the  semi-theoretical  method  we  have  been  using  to  estimate 

dispersion  energies  proved  to  give  values  very  close  to  (variation- 

perturbation )  that  had  been  calculated  using  a  larger  all-electron  basis 
set . 


( 1 ) SDO-MBPT  f  4) 

We  also  calculated  the  EA.  \.Tn  from  correlated  monomer 

tL  ,n 1 P 

multipoles  using  the  larger  all-electron  basis  set.  These  results  suggest 
how  one  might  use  Egj^MTP  and  E^-j^p  results  to  separate  intra  and 

intermolecular  contributions  to  the  correlation  energy. 

V.  Geometry  Optimizations  at  the  MC-SCF/CI  Level 

We  explored  a  number  of  the  features  of  the  GAME5S  program  which  has 
the  capability  of  optimizing  geometries  by  a  derivative  method  at  the  MC-SCF 
and  CASSCF  levels  as  well  as  at  the  SCF  level.  We  are  exploring  and  using 
more  of  the  new  options  in  the  GAMESS  program  for  our  studies  on  the 
energetic  oxetanes. 


02  ' 


VI.  POLY-CRY  ST 


Commensurate  with  the  ONR  priorities  expressed  to  us  by  our  ONR 
Contract  Monitor  to  devote  our  major  attention  to  the  problem  of  the 
cationic  polymerization  of  oxetanes,  we  only  devoted  very  minimal  effort  to 
further  development  and  testing  of  the  POLY-CRYST  program  (for  ab-initio 
calculations  of  crystals  and  polymers). 

For  efficiency  we  initiated  changes  in  the  way  the  integrals  are  picked 
up.  We  also  ran  a  few  test  calculations. 
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I.  Conversion  of  Additional  Programs  to  CRAY  Vector  Supercomputers  and 
Program  Enhancements 

Our  exploratory  ab-initio  MRD-CI  (mul ti reference  double  excitation- 
configuration  interaction)  calculations  on  the  reaction  of  oxetane  with 
protonated  oxetane  for  the  propagation  step  in  cationic  polymerization 
indicated  that  even  for  oxetane  and  certainly  for  the  substituted  oxetanes 
we  were  going  to  have  to  consider  more  active  orbitals,  more  reference 
configurations  and  larger  basis  sets  than  the  MRD-CI  program  was  capable  of 
handl ing  at  that  time . 

Thus,  this  past  year  we  enlarged  these  aspects  of  the  MRD-CI  program 
plus  adding  the  capability  of  taking  more  roots  of  the  Cl  matrix.  The 
programs  are  large  and  complicated  and  such  changes  are  not  trivial.  We 
also  increased  the  capability  for  taking  more  roots  -in  the  explicit  solving 
of  the  Cl  matrix.  This  is  important  if  there  are  potential  energy  surfaces 
involving  excited  states.  We  first  enlarged  these  aspects  of  the  MRD-CI 
program  on  the  NRL  CRAY  XMP-12  which  runs  under  the  COS  operating  system. 
Because  each  such  MRD-CI  calculation  on  the  larger  substituted  energetic 
oxetanes  can  take  significant  amounts  of  computer  time  and  a  number  of 
different  geometries  must  be  calculated  for  each  pair  of  partners,  we  also 
then  added  these  expanded  capabilities  to  the  CTSS  versions  of  the  MRD-CI 
program  which  run  on  the  other  CRAY  computers. 

A  population  analysis  program  to  calculate  charges  on  the  atoms  and 
total  overlap  populations  between  atoms  from  the  density  matrix  following 
the  final  Cl  calculations  was  also  written.  This  will  give  important 
insight  into  such  questions  as  the  relative  charges  on  the  C?  atom  in 

various  substituted  oxetanes  undergoing  reaction  and  the  incipient  making 
and  breaki ng  of  bonds  . 

To  help  in  choosing  molecular  orbitals  most  optimally  localized  in  the 
reaction  (bond  making  and  bond  breaking)  region,  we  wrote  a  program  to 
locate  the  centroids  of  each  localized  molecular  orbital.  We  further 
extended  this  strategy  by  supplementing  the  choice  of  localized  molecular 
orbitals  with  suitable  lone  pairs  or  empty  orbitals  which  could  take  part  in 
making  or  breaking  bonds. 

This  past  year  we  implemented  and  initiated  prototype  calculations  on  a 
major  program  system  GAMESS  (Gaussian  Atomic  and  Molecular  Energy  Systems) 
which  has  the  capability  of  optimizing  geometries  by  a  derivat.ve  method  not 
only  at  the  SCF  level  but  also  at  the  MC-SCF/CI  (multi-configuration 
SCF/conf iguration  interaction)  or  CASSCF  (Complete  Active  Space  MC-SCF) 


levels.  This  can  prove  important  in  any  calculations  on  energetic  compounds 
involving  nitro,  nitrato  and  nitramino  groups.  The  GAMESS  program  was 

2 

originally  written  by  Dr.  Michel  Dupuis  of  the  old  NRC  (National  Resource 
for  Computational  Chemistry).  Dr.  Kaufman  had  been  a  visiting  scientist  at 
2 

NRC  and  had  collaborated  with  Or.  Dupuis  in  using  that  program  for  her 
original  GVB  and  MC-SCF  calculations  on  nitromethane .  With  the 
collaboration  of  Dr.  Victor  Saunders,  the  head  of  the  Atomic  and  Molecular 
Theory  Group  at  Daresbury,  England,  we  got  the  Daresbury  CRAY  COS  version  of 
GAMESS  running  on  the  NRL  CRAY  XMP-12.  Significant  changes  had  been  made  by 
Daresbury  to  the  original  GAMESS  program  to  mesh  in  their  ATMOL  SCF  and 
CASSCF  (Complete  Active  Space  MC-SCF)  with  the  geometry  optimization  scheme. 


As  a  result  of  our  expertise  on  CRAY  computers,  we  acted  as  "friendly 
users"  on  the  NSF  CRAY  XMP-48  supercomputer  at  the  San  Diego  Supercomputer 
Center  (SDSC)  located  at  GTA  San  Diego,  California  for  December  1985  prior 
to  the  official  opening  of  that  center.  That  CRAY  runs  under  a  version  of 
the  CTSS  operating  system. 

As  a  result  of  our  interaction  with  that  SDSC  computer  center,  we  were 
asked  to  serve  on  their  review  committee  for  granting  NSF  computer  time  on 
their  CRAY  XMP-48.  Dr.  P.  C.  Hariharan  served  at  their  January  1986 
computer  allocation  committee  meeting.  Dr.  Walter  S.  Koski  served  at  their 
March,  June,  and  September  meetings  and  Dr.  Joyce  J.  Kaufman  will  serve  at 
their  December  1986  meeting. 


II.  MRD-CI  Calculations  for  the  Propagation  Step  in  Cationic  Polymerization 
of  Energetic  Oxetanes 

Our  major  research  this  past  year  has  been  on  theoretical  quantum 
chemical  ab-initio  MRD-CI  (mul tireference  double  excitation  -  configuration 
interaction)  calculations  on  the  mechanism  of  the  propagation  step  of 
cationic  polymerization  of  oxetane  by  protonated  oxetanes.  Energetic 
polymers  made  by  cationic  polymerization  of  energetic  substituted  oxetanes 
are  of  significant  Navy  interest.  Our  prior  theoretical  quantum  chemical 
calculations  and  the  electrostatic  molecular  potential  contour  (EMPC)  maps 
generated  from  them,  have  enabled  us  to  understand  the  initiation  step  and 
to  predict  the  propensity  of  different  energetic  substituted  oxetanes  to 
undergo  cationic  polymeri zation  and  to  pick  optimal  copolymer  partners  even 
prior  to  synthesis  of  the  monomers  themselves.  Our  current  MRD-CI 
calculations  will  shed  insight  on  the  subsequent  propagation  step  in 
cationic  polymerization  and  how  different  copolymer  partners  would  prefer  to 
propagate.  This  influences  the  subsequent  final  polymer. 

The  oxetane  -  protonated  oxetane  systems  (especially  the  substituted 
ones)  are  larger  than  any  molecular  reacting  systems  for  which  any  MRD-CI 
calculations  have  been  carried  out.  We  developed  and  implemented  a  new 
approach  based  on  localized  occupied  and  virtual  orbitals  in  the  interaction 
region  with  the  remainder  of  the  localized  molecular  orbitals  being  folded 
into  an  effective  Cl  Hamiltonian.  This  method  is  completely  general  and 
applicable  to  reactions  and  molecular  decompositions  of  energetic  compounds 
for  ground  and  electronically  excited  states. 

After  initiation  the  succeeding  step  in  cationic  polymerization  is 
propagation.  In  the  past  some  experimentalists  had  suggested  that  there 
would  be  ring  opening  of  the  protonated  carbocation 


H 


which  would  then  be  the  propagating  species.  The  experimental ists  had 
further  suggested  that  the  more  electrophilic  a  protonated  carbocation  was, 
the  more  efficient  a  propagating  species  it  would  be.  Our  previous  quantum 
chemical  calculations  had  shown  that  by  any  criterion,  the  protonated 
carbocation  of  3,3-di nitrooxetane  would  be  the  best  propagating  species. 
Conversely,  our  EMPC  maps  had  shown  that  3 ,3-di ni trooxetane  would  have  the 
least  propensity  to  undergo  initiation  of  cationic  polymerization.  Since 
3 , 3-di n i trooxetane  does  net  undergo  cationic  pol ymer i zat i on  under  any  usual 
conditions,  it  is  the  initiation  step  which  is  dominant  in  the  propensity 
for  cationic  po 1 ymer i zat i on  . 
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However,  experimentalists  believe  that  1°  carbocations  do  not  form 
except  under  very  rare  circumstances  and  are  not  a  very  stable  species  and 
thus  would  not  be  expected  to  exist  long  enough  in  aqueous  solution  to  serve 
as  the  propagating  species. 

Recently  (June  1985)  several  different  experimentalists  in  cationic 
polymerization  suggested  to  us  that  a  possible  mechanism  would  be  attack  of 
protonated  oxetanes  on  oxetanes  (or  vice  versa)  with  concomitant  ring 
opening  of  the  protonated  oxetane  according  to  the  following  general  scheme 


H 


Manzer  suggested: 


lillya  suggested  essentially  the  same  overall  picture: 


Thus  we  initiated  a  major  research  effort  on  ab-initio  MRD-CI 
calculations  on  opening  the  ring  of  oxetane  or  protonated  oxetane  and  on  the 
opening  of  the  protonated  oxetane  ring  in  the  course  of  interaction  with 
oxetane.  MRD-CI  calculations  are  necessary  if  a  bond  is  being  broken  or 
formed  since  often  there  are  several  Slater  determinants  involved  in  the 
electronic  wave  functions.  We  have  already  carried  out  the  ab-initio  MRD-CI 
calculations  for  a  number  of  different  C  -  0  interring  distances  and  angles 
and  different  angles  of  opening  of  the  oxetane  and  protonated  oxetane  ring. 
The  substituted  oxetanes  are  large  molecular  systems  and  their  interactions 
with  protonated  substituted  oxetanes  lead  to  even  larger  systems.  These  are 
larger  molecular  systems  than  have  ever  been  calculated  with  MRD-CI  methods 
and  also  exceed  the  amount  of  data  that  can  be  handled  even  in  the  current 
CRAY-XMP  series.  Hence,  we  developed  and  implemented  a  new  MRD-CI  approach 
based  on  localized  orbitals  occupied  and  virtual  in  the  reaction/interaction 
regaion  with  the  remainder  of  the  nonparticipating  localized  occupied 
molecular  orbitals  being  folded  into  an  effective  Cl  Hamiltonian. 

A.  Validation  of  Localized  Orbital  MRD-CI  Technique 

Many  years  ago,  we  derived  and  implemented  a  technique  in  which  the 
effects  of  the  occupied  molecular  orbitals  from  which  excitations  are  not 

allowed  are  folded  into  an  effective  Cl  Haniltonian.1  This  has  the 
advantage  that  it  is  no  longer  necessary  to  transform  the  integrals  over 
atomic  orbitals  to  integrals  over  molecular  orbitals  from  which  excitations 
are  not  allowed.  It  is  this  transformation  which  is  the  computer  resource 
(space-  and  time-wise)  intensive  part  of  the  MC-SCF  or  especially  of  the  Cl 
probl em. 

Use  of  the  effective  Cl  Hamiltonian,  employing  sets  of  highest  occupied 
(and  lowest  unoccupied)  conventional  delocalized  molecular  orbitals,  is 
applicable  to  and  has  been  used  previously  for  calculations  of  such  physical 
phenomena  as  ionization  potentials,  electron  affinities  and  electronic 
spectra.  However,  for  molecular  decompositions  and  molecular  reactions,  all 
of  the  molecular  orbitals  in  the  spatial  region  of  interest  will  not  be 
among  the  highest  occupied  and  lowest  unoccupied  orbitals.  Thus,  we 
2 

proposed  in  1980,  for  molecular  decompositions  and  reactions  to  transform 
to  localized  molecular  orbitals  and  to  carry  out  the  Cl  calculations 
explicitly  on  the  localized  molecular  orbitals  occupied  and  virtual  in  the 
spatial  region  of  interest  folding  the  remainder  of  the  occupied  localized 

2  3 

molecular  orbitals  into  an  effective  Cl  Hamiltonian.  ’ 

We  showed  by  test  examples  (with  a  system  of  such  a  size  that  the 
entire  valence  space  MRD-CI  could  be  run  on  the  NRL  CRAY-XMP)  that  the  MRD- 
CI  based  on  localized  orbitals  gave  a  potential  energy  surface  for  molecular 
decomposition  essentially  parallel  to  that  using  the  entire  valence  space 
MRD-CI.  The  test  example  that  we  chose  was  the  H^C-F  dissociation  of 

f 1 uoromothane  using  the  full  valence  space  and  using  only  localized  orbitals 
in  the  dissociation  region  (folding  the  remainder  of  the  occupied  molecular 
orbitals  into  an  effective  Cl  Hamiltonian).  These  calculations  were  carried 
out  for  several  purposes: 
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1.  To  validate  our  hypothesis  that  potential  energy  surfaces  for  bond 
dissociations  from  MRD-CI  calculations  based  on  strictly  orthogonal 
localized  molecular  orbitals  in  the  bond  dissociation  region  (with  the 
remainder  of  the  occupied  orbitals  folded  into  an  effective  Cl  Hamiltonian) 
would  be  essentially  parallel  to  MRD-CI  potential  energy  surfaces  calculated 
using  full  valence  molecular  orbital  space.  For  this  purpose  we  used  a 

standard  Dunning  contraction  of  the  Huzinaga  basis  set  (9s5p->  3S2P)  on  C  and 

F  and  a  3s  -  Is  basis  set  on  the  hydrogens  and  ran  MRD-CI  calculations 
allowing  excitations  from  all  the  occupied  valence  molecular  orbitals. 

With  this  same  basis  set,  we  then  carried  out  the  MRD-CI  calculations 
using  the  localized  molecular  orbitals  in  the  C  -  F  bond  and  on  the  F  plus 
the  full  valence  virtual  localized  mol°cuiar  orbital  space. 

2.  To  show  that  MRD-CI  potential  energy  surfaces  calculated  with  the 
all-electron  basis  set  with  the  same  exponents  and  contraction  coefficients 
as  our  customary  ab-initio  M0DP0T  valence-electron  basis  set  (incorporating 
ab-initio  effective  core  model  potentials  for  inner  shells)  were  essentially 

parallel  to  those  calculated  with  a  more  standard  Dunning  9S5P+  3S2P 
contraction . 

The  atomic  all-electron  basis  set  for  this  set  of  MRD-CI  calculations 

was  10s6p-»  2S1P  on  C  and  F  and  the  same  3s->  Is  atomic  basis  set  on  H  as 
above.  The  SCF  calculations  were  carried  out  and  the  canonical  molecular 
orbitals  were  transformed  to  strictly  orthogonal  localized  molecular 
orbitals.  The  MRD-CI  calculations  were  carried  out  both  with  the  full 
valence  molecular  orbital  space  (occupied  and  virtual)  and  then  with  the 
localized  molecular  orbital  in  the  C  -  F  bond  plus  the  full  valence  virtual 
orbital  space. 

The  four  different  potential  energy  surfaces  are  shown  in  Figure  1. 


Fluoromethdne  MRD-CI  Energy  (Extrapolated) 


x 


full  spoct 
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From  Figure  1  it  can  be  seen  that: 

1.  The  MRD-CI  potential  energy  surfaces  are  essentially  parallel 

for  the  9S5^->  3S2^  basis  set  both  for  the  full  all-valence  electron  space 
and  for  the  restricted  localized  orbital  space. 

2.  The  MRD-CI  potential  energy  surfaces  are  essentially  parallel 

for  the  all-electron  MODPOT  1 0s 6^-  2S1^  basis  set  both  for  the  full  all 
valence  electron  space  and  for  the  localized  orbital  space. 

3.  The  MRD-CI  curves  are  also  essentially  parallel  for  the 
Dunning  3s 2^  and  the  all-electron  MODPOT  2S1^  basis  set. 

Thus,  the  localized  orbital  MRD-CI  approach  appears  valid  and  well 
behaved.  This  technique  holds  great  promise  for  calculations  on  larger 
systems . 


B.  MRD-CI  Cal cu1 ations  on  the  Reaction  of  Oxetane  with  Protonated  Oxetane. 


These  MRD-CI  calculations  for  the  reaction  of  protonated  oxetane  with 
oxetane  are  a  computational 1y  and  labor  intensive  project.  For  each 
different  inter-  and  intra-molecular  geometry  point,  first  the  SCF 
calculation  must  be  run,  then  the  resulting  SCF  canonical  delocalized 
molecular  orbitals  must  be  localized.  We  wrote  a  program  to  determine  the 
centroids  of  the  localized  bonds  compared  to  the  bonds  and  atoms  involved  in 
the  interaction/reaction.  We  found  we  had  to  supplement  that  choice  with 
the  empty  orbitals  and  lone  pairs  that  might  be  involved  in  the  bond  making 
and  bond  breaking. 


Next  a  small  single  reference  Cl  calculation  must  be  carried  out  to 
determine  the  major  reference  configurations.  Then  a  larger  MRD-CI 
calculation  is  run.  The  configurations  corresponding  to  all  of  the  single 
and  double  excitations  are  then  generated  relative  to  the  multireference 
configuration  wave  function.  The  program  estimates  the  energy  contribution 
of  each  of  these  single  and  double  excitation  configurations  by  a 
perturbation  procedure.  All  of  the  configurations  contributing  more  energy 
than  a  predetermined  threshold  are  included  explicitly  in  the  Cl  calculation 
(E^).  Then  the  energies  of  the  excluded  other  configurations  are  added 


back  in  by  an  extrapolation  perturbation  procedure 


(eext 


)  and  final 


1  y  a 


Davidson  type  correction  is  added  to  correct  for  size  extensivity. 


ref  2 

E ( f u 1 1  Cl  estimate)  =  E  (EXT)  +  (1-  £  c  )  [E(EXT)  -  E(Ref)] 

P  P 

The  MRD-CI  method  used  is  that  of  Buenker  and  Peyerimhoff  and  the 

definitions  of  the  various  energies  are  consistent  with  those  used  by 

..  4-7 

them. 


We  reported  last  year  that  we  had  carried  out  test  MRD  Cl  calculations 
on  oxetane  and  protonated  oxetane  and  for  their  ring  opened  species  using  a 

standard  Dunning  9s5p-»  3S2P  all-elcctron  basis  set  and  our  ab-initio  MODPOT 
basis  set.  The  MRD-CI  results  for  energy  differences  showed  the  same  trends 
for  both  basis  sets.  Moreover,  the  difference  in  the  ab-initio  MODPOT  MRD- 
CI  energies  between  (at  all  oxetane  and  protonated  oxetane  (at  all  Cl 
levels)  levels  E^,  E^y,  was  a-u-  (8-43  ev.)  in  excellent 

Q 

agreement  with  the  experimental  proton  affinity  of  oxetane  (8.56  ev.)  . 

The  proton  affinity  contributes  sufficient  energy  to  influence  the 
protonated  oxetane  ring  to  open. 

Thus,  all  of  the  calculations  on  the  oxetane  plus  protonated  oxetane 

[which  we  call (double  oxetane)+]  were  carried  out  with  our  ab-initio 
MODPOT/VRDDO  technique  meshed  in  with  the  MRD-CI  program. 

For  the  MRD-CI  calculations  on  oxetane  and  protonated  oxetanes 
themselves  (and  their  ring  opened  species)  (Figure  2), 
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we  considered  the  localized  orbitals  in  the  C^-O^  bond,  the  C-^-C^  bond,  the 

C^-C^  bond,  the  Oj-H  +  bond,  the  lone  pairs  in  Oj  and  the  bonds  connecting 
hydrogens  to  C2  and  C^. 


We  first  ran  a  modest  size  Cl  calculation, 
contributions  of  the  various  configurations  to 
point.  From  these  Cl  results  we  picked  the  mos 
configurations,  generated  all  single  and  double 
all  the  reference  configurations,  estimated  by 
energy  contributions  of  each  configuration,  set 
carried  out  the  MRD-CI  calculations,  added  back 
contributions  of  configurations  not  included  ex 
extrapolated  to  get  and  then  added  the  Dav 

account  for  size  extensivity. 


We  monitored  the 
the  wave  function  at  each 
t  important  reference 
excitations  with  respect  to 
a  perturbative  procedure  the 
a  threshold  (0.000035  a.u.), 
in  the  perturbative  energy 
pi  icitly  in  the  MRD-CI , 
idson  type  correction  to  help 


There  were  several  significant  results  on  the  single  ring  compounds: 
1.  While  for  the  closed  ring  oxetane  the  ground  state  SCF  wave 

function  contributed  -96%  of  the  final  MRD-CI  wave  function,  for 


the  ring  opened  oxetane  (  S  -  19u)  the  ground  state  SCF  wave 
function  contributed  only  70%  of  the  final  MRD-CI  wave  function. 
Conversely,  however,  for  the  closed  ring  protonated  oxetane,  the 
SCF  wave  function  contributed  -  97%  °f  the  final  MRD-CI  wave 
function.  Conversely,  however,  for  the  closed  ring  protonated 
oxetane  the  SCF  wave  function  contributed  -97%  of  the  final  Cl  MRD- 
CI  wave  function.  Analysis  of  the  final  wave  function  indicates 


-v\>VA 
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that  in  the  case  of  the  protonated  oxetane  both  electrons  remain 
with  the  oxygen  atom  when  the  C^-0  bond  opens. 

2.  Inspection  of  the  gross  atomic  populations  (both  SCF  and  MRD-CI)  of 
protonated  oxetane  (Table  1)  indicates  that  there  is  not  a  positive 
charge  on  the  protonated  oxygen  but  a  negative  charge.  This  is  due 
to  the  electronegativity  of  oxygen  being  larger  than  that  of  any  of 
the  other  atoms.  There  are  also  varying  degrees  of  negative  crarge 
on  the  carbon  atoms  in  protonated  oxetane  since  the  carbons  are 
more  electronegative  than  the  hydrogens.  The  positive  charges 
reside  on  the  hydrogen  atoms.  In  oxetane  itself  the  oxygen  atom  is 
somewhat  less  negatively  charged  than  in  protonated  oxetane  and  the 
carbons  have  varying  degrees  of  negative  charge  while  there  are 
positive  charges  on  the  hydrogen  atoms.  This  charge  redistribution 
has  a  profound  influence  on  how  one  should  consider  the  mechanism 
of  the  reaction  between  oxetanes  and  protonated  oxetanes  in  the 
propagation  step  of  cationic  polymerization.  In  Table  I  the  GAP's 
for  the  ring-opened  oxetane  and  ring-opened  protonated  oxetane  are 
presented  also  at  the  SCF  and  MRD-CI  levels.  The  MRD-CI  GAP  on  the 
oxygen  of  ring-opened  oxetane  is  somewhat  smaller  than  the  GAP  from 
the  SCF  result.  This  is  due  to  the  result  that  the  ground  state 
SCF  wave  function  of  ring-opened  protonated  oxetane  only 
contributes  70%  to  the  final  MRD-CI  wave  function. 


17 


We  next  carried  out  similar  MRD-CI  calculations  for  the  cationic 
polymerization  ring  opening  of  protonated  oxetanes  by  oxetane  (Figure  3). 
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We  have  carried  out  calculations  for  a  great  many  geometries  for 
various  paths  of  attack  for  the  system  oxetane  plus  protonated  oxetane.  We 
have  varied  the  angle  of  opening  of  the  protonated  oxetane  ring  (ring  'A'), 
varied  the  internuclear  distance  between  the  rings  and  angle  between  the 
plane  of  the  oxetane  ring  (ring  '3')  and  the  plane  of  the  protonated  oxetane 
ring. 
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GEOMETRY  VARIATIONS  FOR  OXETANE  (RING  B) 
PLUS  PROTONATED  OXETANE  (RING  A) 


A  'A'  ring  6=0°  (Closed  Ring) 

B  'A*  ring  o  =  5° 

E  'A'  ring  «  19l  (Open  Ring) 

F  'A'  ring  like  propanol 

G  'A'  ring  propanol  with  H5A,  H6A  dihedral  angles  rotated  by  180°  to 

i  70°  (away  from  'A1  ring) 

H  Same  as  1 G '  but  H5A,  H6A  ideally  tetrahedral 

I  Likp  1 E  *  6  »  19°  H5A,  H6A  defined  as  in  H 

GEOMETRY  -  INTERNUCLEAR  BETWEEN  RINGS 
2.4,  2.9,  3.4,  3.9,  4.4,  4.9  a.u. 


a  value  dihedral  angle  between  rings  90°,  135°,  180°  (coplanar) 

Y  109.5°  tetrahedral  around  01G 
132.94°  planar  around  01B 

0  angle  of  bonds  around  C4A 
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If  oxetane  approaches  protonated  oxetane  as  shown  in  Figure  3,  our  MRD- 
CI  calculations  indicate  the  interaction  would  be  repulsive  because  of  the 


interference  of  the  2  H  atoms  on  C, 


Thus  we  have  also  carried  out  MRD-CI 


calculations  pulling  the  positions  of  these  two  H  atoms  back  from  the 
approach  of  the  attacking  oxetane.  as  indicated  below  (Figure  4) 


In  discussions  with  several  JHU  organic  chemists  who  are  experts  in  the 
field  of  carbonium  (carbocation )  chemistry  it  was  further  suggested  to  us 


that  a  possible  geometric  pathway  would  be  similar  to  an  SN2  attack.  An 


attack  would  have  a  similar  pattern  of  five  bonds  (or  incipient  bonds)  to  a 
carbon  in  the  intermediate.  This  would  imply  a  general  mechanism  (Figure  5) 


I  *  „H 

X  — x  *  x — c  — 


C —  H 
"*  H 


(where  X  is  a  closed  shell  species  with  an  electron  pair  pointing  toward 
the  C)  with  an  inversion  (a  Walden-type)  inversion).  In  such  a  case  an 


attack  of  X  would  be  most  favored  if  the  attack  of  X  were  linear  along 
the  C  -  X  bond  direction.  This  mechanism  applied  to  the  oxetane  plus 
protonated  oxetane  system  with  concomitant  opening  of  the  protonated  oxetane 
ring  would  suggest  that  a  possible  favorable  attack  would  be  for  the  0,„  of 


ring  would  suggest  that  a  possible  favorable  attack  would  be  for  the  0^B  of 
the  oxetane  ring  to  attack  protonated  oxetane  along  the  line  of  the  C^-O^ 
bond.  We  are  carrying  out  MRD-CI  calculations  for  this  inode  of  attack.  Our 


preliminary  MRD-CI  results  appear  favorable  for  forming  a  bond  between  Ojg 
and  C^with  concommitant  opening  of  the  C^-0^  bond.  (Figure  6) 
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Examination  of  the  charge  distribution  of  the  oxetane  plus  protonated 
oxetane  systems  shows  that  in  protonated  oxetane  there  is  not  a  positive 
charge  on  the  0  atom  -  but  rather  a  negative  charge.  Also  there  are 
negative  charges  on  the  ring  carbon  atoms.  The  positive  charges  reside  on 
the  hydrogen  atoms  (both  the  H  atoms  on  ring  carbons  and  the  H  on  the  0). 
In  the  intermediate  or  product 
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there  are  still  no  positive  charges  on  either  of  the  oxygens  or  on  the 
carbons  (Table  II)  This  is  somewhat  different  than  the  charge  distribution 
picture  depicted  to  us  by  the  experimentalists.  This  could  have  significant 
implications  on  visualizing  the  mechanism  of  the  propagation  step. 


Table  II 


OXETANE  PLUS  PROTONATEL)  OXETANE 
Ab-Initio  M0D0T/VRDD0 

Gross  Atomic  Populations  -  Valence  Electrons  Only 
RQr  -  2.9  bohr  6  =  19'  y  =  132.94  01A-C4A-01 3  Coll  inear 


SCF 


MRD-CI 


The  total  overlap  populations  (TOP'S)  of  the  oxetane  plus  protonated 
oxetane  system  for  a  particular  path  of  attack  (Table  III)  clearly  indicate 
that  for  the  system  where  the  protonated  oxetane  ring  ('A'  ring)  is  open 

(6=19°)  that  the  inter-ring  b-C4A  9ets  ^ar9er  ^rom  an  ^1B~^4A  inter” 


ring  distance  of  4.4  bohrs  up  until  the  energy  minimum  at  2.9  bohrs  and  then 
starts  to  get  smaller.  This  indicates  formation  of  an  0^-C^  bond. 


Conversely,  when  the  'A'  ring  is  closed  (6*0°)  the  TOPgi b-C4A  more  and 

more  strongly  negative  (indicating  repulsion)  as  the  rings  come  closer.  The 
intra-ring  T0pgiA-C4A  of  the  closed  protonated  oxetane  ring  becomes  less  and 

less  bonding  as  the  oxetane  approaches  indicating  that  the  protonated 
oxetane  ring  will  tend  to  open  upon  attack  by  the  oxetane.  For  many  years 
we  have  found  TOP'S  to  be  an  excellent  sensitive  indicator  of  mechanisms  of 
reactions . 
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Table  III 


OXETANE  PLUS  PROTONATEO  OXETANE 
TOTAL  OVERLAP  POPULATIONS 
VALENCE  ELECTRONS  ONLY 


'A'  Ring 
6  -  0° 

Fully  Closed 
TOP 


'A'  Ring 
6  -  19° 
Fully  Open 
TOP 


R01B-C4A 

(BOHRS) 

01B-C4A 

01A-C4A 

01B-C4A 

01A-C4A 

2.4 

-0.49 

-0.05 

0.39 

-0.03 

2.9 

-0.44 

0.15 

0.42 

-0.02 

3.4 

-0.28 

0.28 

0.37 

-0.01 

3.9 

-0.13 

0.36 

0.27 

-0.01 

4.4 

-0.06 

0.39 

0.17 

-0.00 

a  »  180‘ 


132.94 


0 1 A-C4A-0 1 o  COLLINEAR 


We  have  carried  out  calculations  for  a  variety  of  different  geometries 
as  indicated  above.  The  most  favorable  pathway  we  have  found  for  oxetane 
attacking  protonated  oxetane  is  the  one  suggested  to  us  by  the  JHU  organic 
experimentalists,  namely,  that  the  O^g  (of  oxetane)  will  prefer  to  attack 

linearly  along  the  bond  (of  protonated  oxetane)  with  concomitant 

opening  of  the  bond. 

For  the  MRO-CI  calculation  at  each  point,  we  monitor  the  contributions 
of  the  dominant  and  other  possibly  significant  determinants  to  the  MRD-CI 
wave  functions. 

Some  representati ve  results  of  the  MRD-CI  calculations  are  shown  in 
Tables  IV-XV.  In  these  tables  we  list  the  number  of  configuration  state 
functions  generated  and  the  number  of  these  used  in  each  step,  the  energies 

?  o 

(Escf»  Eqj,  ar,d  E pyLL )  *  the  ^  sur;i  c<"  a1'*  tlie 

2 

reference  congf igurations )  and  the  c  for  each  function  contributing  more 

p 

than  c  "=  0.005  to  the  final  MRD-CI  wave  function.  We  are  now  initiating 
similar  studies  on  substituted  oxetanes  and  protonated  oxetanes  to 
understand  their  preferences  for  reacting  with  themselves  or  with  other 
substituted  species.  This  will  give  insight  into  copolymer  propensities. 
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OXETANE  PLUS  PROTONATED  OXETANE 


Ab-Initio  MODPOT/VRDDO  MRD-CI 

'A'  ring  6-0°  r  =  2.4  bohrs  a  -  180°  y  -  182.94°  01A-C4A-01B 

fully  closed  collinear 


ENERGIES  (a.u.) 

escf 

-70.751664 

ECI 

-70.945929 

(891  csf's  at  10  MH) 

eext 

-70.953768 

(30157  csf's  generated) 

efull 

-70.957530 

Ic2 

0.9617 

GROUND  STATE 

0.9070 

(01A-H0+)2  — > 

(01A-H0+*)2 

0.0065 

(01B-C2B) 2  —  > 

(01B-C2B*)2 

0.0086 

(C2A-01A) 2  — > 

(C2A-01A*) 2 

0.0097 

(01B-C4B) 2  — > 

(01B-C4S*) 2 

0.0073 

(C3A-C4A)2  — > 

(C3A-C4A*)2 

0.0064 

(C4A-H5 A ) 2  — > 

(C4A-H5A*)2 

0.0081 

(C4A-H6A)2  - > 

(C4A-H6A*) 2 

0.0081 

K.‘- 
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OXETANE  PLUS  PROTONATED  OXETAHE 
Ab-Initio  MODPOT/VRDDO  MRD-CI 


'A'  ring  5  »  0°  r  =  2.9  bohrs  a  -  180°  y  =  132.94- 
fully  closed 


-SCF 


‘Cl 


-EXT 


FULL 


GROUND  STATE 

,2 


(01A-H0+)  — >  (01A-H0+ 


(C2A-01A) c  — >  (C2A-01A 


*,  2 


(01 8-C4B) ^  — >  (01B-C4C 


*  7 


(C4A-01Ap  — >  (C4A-01A  )4 


(01B-C2BP  — >  (01B-C2B 


(C4A-H6A 


1 2 


— > 


14A-H6A*)2 


(C3A-C4A )  ‘ 


->  (C3A-C4A 


*,  2 


(C4A-H5AP  — ->  (C4A-H5A 


01A-C4A-01B 
col  1  inear 


ENERGIES  (a.u.) 
-71.077815 


-71.267822 

(794  csf's  at  10  pH) 

-71.274607 

(33881  csf's  generated) 
-71.277755 


0.9639 


0.9074 

0.0060 

0.0085 

0.0091 

0.0021 

0.0101 

0.0071 

0.0067 

0.0071 
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Table  VI 


OXF.TANE  PLUS  PROTONATEU  OXETANE 
Ab-Initio  MODPOT/VRDDO  MRD-C J 


'A'  ring  6-0°  r  -  3.4  bohrs 
fully  closed 


130°  f  -  132.94' 


GROUND  STATE 


01A-C4A-01 B 
col  1  inear 


ENERGIES  (a.u.) 
-71.230575 

-71.415443 

(692  csf's  at  10  pH) 

-71.421751 

(33881  csf's  generated) 
-71.424386 

0.9665 

0.9073 


(01A-H0+)2 

— ->  (01A-HO+*)2 

0.0057 

(C3A-C4A)2 

— >  (C3A-C4A*)- 

0.0068 

(C2A-01A)2 

- >  (C2A-01A*) 2 

0.0078 

(C4A-01A) 2 

- >  (C4A-01 A*) 2 

0.0039 

(01B-C4B)2 

- >  (01B-C4B*)2 

0.0105 

(01B-C2B)2 

- >  (01B-C2C*)2 

0.0110 

(C4A-H5A) 2 

— >  (C44-H5A*) 2 

0.0068 

(C4A-H6A) 2 

- >  (C4A-H6A*) 2 

0.0068 

’VV’m*  *«,►  V  V  *«*  V  V  *■* "\f"V  J*  'V1, 


table  VII 


OXETANE  PLUS  PROTONATED  OXETANE 
Ab-Initio  MOOPOT/VRDDO  MRD-CI 


'A'  ring  6  »  0°  r  -  4.4  bohrs  a 
fully  closed 

ESCF 

ECI 

eext 

E  FULL 

Ic2 

GROUND  STATE 
(01A-H0+)2  --> 

(C4A-H6A ) 2  ---> 

(C4A-01 A  )  2  — -> 

(C4A-H5A)2  ---> 

(C3A-C4A) 2  — -> 

(01B-C2B)2  — -> 

(01 B-C4B) 2  — -> 


180°  y  =  132.94°  01A-C4A-01 B 

col  1  inear 

ENERGIES  (a.u.) 
-71.313176 

-71.492765 

(545  csf's  at  10  pH j 

-71.498149 

(33881  csf's  generated) 
-71.500286 

0.9696 

0.9072 

0.0054 

0 . 0066 

0.0059 

0.0066 

0.0070 

0.0120 


(01fl-H0  +  V 

(C4A-H6A*) 2 
(C4A-01A*)2 
(C4A-H5A*) ? 
(C3A-C4A*)2 
(01B-C2B*)2 
(0) B-C4B* j2 


0.0119 


$9 
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Table  VIII 


OXETANE  PLUS  PROTONATED  OXETANE 
Ab-Initio  MODPOT/VRDDO  MRD-CI 


'A'  ring  6-19°  r  -  2.4  bohrs  c  -  180°  y  -  132.94 
fully  open 


0  01A-C4A-01B 

col  1  inear 


GROUND  STATE 

(01A-HO+)2  — >  (01A-H0+*) 2 

(01B-C4B)2  — >  (01B-C4B*)2 

(C3A-C4A)2  — ->  (C3A-C4A*)2 

( 0 1 B-C4A ) 2  --->  (C4A-01B*)2 

(01B-C2B)2  — ->  (01 B-C  2  3* ) 2 

(C2A-01 A) 2  — >  (C2A-01A*) 2 

(C4A-H6A)2  — ->  (C4A-H6A*)2 

(C4A-H5A) 2  — >  (C4A-H5A*) 2 


ENERGIES  (a.u.) 
-71.267994 


-71.436155 

(546  csf's  at  10  pH) 

-71.440836 

(33881  csf's  generated) 
-71.442487 


0.9742 


0.9164 


0.0080 


0.0069 


0.0066 


0.0028 


0.0071 


0.0125 


0.0068 


0.0069 


V*  v*  wuvuv 
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OXETANE  PLUS  PROTONATED  OXETANE 
Ab-Initio  MODPOT/VRDDO  MRD-CI 


'A'  ring 

6  -  19°  r  -  2.9  bohrs  a  -  180° 

y  -  132.94°  01A-C4A-01B 

fully  open 

col  1  inear 

ENERGIES  (a.u.) 

escf 

-71.342800 

ECI 

-71.515734 

(520  csf's  at  10  pH) 

EEXT 

-71.520481 

(33881  csf's  generated) 

efull 

-71 . 522302 

Ic2 

0.9720 

GROUND 

STATE 

0.9112 

(01A-H0+)2 

— ->  (01A-HO+V 

0.0078 

(C4A-H5A)2 

— >  (C4A-H5A*)2 

0.0067 

(01B-C4B)2 

- >  (01 B-C4B*) 2 

0.0073 

(01 8-C  2  3 ) 2 

- >  (01 B-C2B*) 2 

0.0074 

(C2A-01A) 2 

— >  (C2A-01A*)2 

0.0124 

(01C-C4A)2 

- >  (C4A-01B*)2 

0.0058 

(C3A-C4A) 2 

— ->  (C3A-C4A*) 2 

0.0067 

(C4A-H6A) 2 

- >  (C4A-H6A*) 2 

0.0067 

a  •  vT ^  i  - 
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Table  X 
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OXF.TANE  PLUS  PROTONATED  OXETANE 
Ab-Initio  MODPOT/VRDDO  MRD-C I 


'A'  ring  <5-19°  r  -  3.4  bohrs  a  -  180°  y  -  132.94°  01A-C4A-01B 

fully  open  collinear 


"SCF 


’Cl 


'EXT 


'FULL 


ENERGIES  (a.u.) 
-71.322213 

-71.497370 

(510  csf's  at  10  pH) 

-71.501523 

(33881  csf's  generated) 
-71.503595 


Ic2 

GROUND  STATE 

0.9689 

0.9058 

(01A-H0+)2  — > 

(01A-H0+*)2 

0.0077 

(C2A-01A) 2  — > 

(C2A-01 o*) 2 

0.0122 

(C4A-H5A)2  — > 

(C4A-H5A*) 2 

0.0066 

(02B-C2B)2  — > 

(02B-C2B*) 2 

0.0079 

(01B-C4I3) 2  — > 

(01B-C4B*)2 

0.0078 

(C4A-H6A) 2  — > 

(C4A-H6A*)2 

0.0066 

(C3A-C4A)2  — > 

(C3A-C4A*) 2 

0.0066 

(01B-C4A)2  — > 

(C4A-01B*)2 

0.0077 

Table  XI 


OXETANE  PLUS  PROTONATED  OXETANE 
Ab-Initio  MODPOT/VRDDO  MRD-CI 


'A'  ring  6-19°  r  -  4.4  bohrs  a  =  180°  y  -  132.94° 
fully  open 


E 


SCF 


ENERGIES  ( 
-71.250019 


ECI  -71.419391 

(470  csf's 

Eext  -71.423061 

(33881  csf 

EpULL  -71.425320 


lcZ  0.9650 

GROUND  STATE  0.9037 


(01 B-C4B) 2 

- > 

(01B-C4B*)2 

0.0092 

(01A-HQ+)2 

- > 

(01A-H0+*) 2 

0.0075 

(C4A-H6A) 2 

- > 

(C4A-H6A*)2 

0.0065 

(C4A-H5A)2 

- > 

(C4A-H5A*)*' 

0.0065 

(C2A-01A)2 

- > 

(C2A-01A*)2 

0.0121 

(C3A-C4A) 2 

- > 

(C3A-C4A*) 2 

0.0064 

(01B-C2B)2 

- > 

(018-C2B*)2 

0.0093 

(01 B-C4A ) 2 

- > 

(C4A-01B*)2 

0.0036 

I 

01A-C4A-01B 
col  1  inear 

at  10  pH) 

s  generated) 


Table  XII 


OXETANE  PLUS  PROTONATED  OXETANE 
Ab-Initio  MODPOT/VRDDO  MRD-CI 

*A'  ring  6-19°  r  -  2.4  bohrs  a  =  135°  y  -  132.94°  01A-C4A-01B 

fully  open  collinear 

ENERGIES  (a.u.) 

Escf  -71.269251 


EEXT 

efull 


-71.437514 

(562  csf's  at  10  pH) 

-71.441919 

(33881  csf's  generated) 
-71.443578 


lc2  0.97^1 

GROUND  STATE  0.9164 


(01A-H0+ ) 2 

- > 

(01A-H0+V 

0.0080 

(C4A-01B)2 

- > 

(C4A-0I3*)2 

0.0028 

(01B-C2B)2 

- > 

(01B-C2B*)’ 

0.0071 

(01B-C4B)2 

—  > 

(01B-C4B*)2 

0.0070 

(C4A-H5A)2 

- > 

(C4A-H5A*) 2 

0.0067 

(C2A-01 A) 2 

- > 

(C2A-01A*)2 

0.0125 

(C3A-C4A)2 

- > 

(C3A-C4A*)2 

0.0068 

(C4A-H6A ) 2 

- > 

( C4A-H6A* ) 2 

0.0069 

V.'. 

W 

A"- 
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OXF.TANE  PLUS  PROTONATEO  OXETANE 
Ab-Initio  MODPOT/VRDDO  MRD-CI 


'A'  ring  6  -  19°  r  =  2.9  bohrs  u  =  135J  y  =  132.94“ 
fully  open 

ENERGIES  (a.u 

Escp  -71.343003 


EEXT 

efull 


-71.515998 
(528  csf's  at 

-71.515998 
(33881  csf's 

-71.522023 


Ic2  0.9721 

GROUND  STATE  0.9112 


(01A-H0+) 2 

- > 

(01A-H0+*)2 

0.0078 

(C4A-H5A)2 

- > 

(C4A-H5A*) 2 

0.0066 

(01 3-C43) 2 

- > 

(013-C4B*)2 

0.0073 

(013-C2B)2 

- > 

(01 B-C2B*) d 

0.0074 

(C2A-01A)2 

- > 

(C2A-01A*)2 

0.0124 

(C4A-H6A)2 

- > 

(C4A-HSA*) 2 

0.0067 

(01 3-C4A) 2 

- > 

(C4A-01B*)2 

0.0058 

(C3A-C4A)2 

- > 

(C3A-C4A*) 2 

0.0067 

01A-C4A-01 B 
col  1  inear 

J. 

10  pH) 
generated) 


Table  XIV 


OXETANE  PLUS  PROTONATED  OXETANE 
Ab-Initio  MODPOT/VRDDO  MRD-CI 


’A'  ring  6-19°  r  -  4.4  bohrs  u  =  135°  y  “  132.94°  01A-C4A-01B 

f ul ly  open  col  1  inear 


ENERGIES  (a.u, 
-71.249965 


-71.419232 

(463  csf's  at  10  pH) 

-71.423010 

(33881  csf's  generated) 


efull 

-71.425265 

Ic2 

0.9651 

GR0UN0  STATE 

0.9038 

(01B-C4B) 2 

- > 

(01B-C4B*)2 

0.0093 

(01A-HO+)2 

- > 

(C1A-H0+*) 2 

0.0075 

(C4A-H5A)2 

- > 

(C4A-H5A*)2 

0.0065 

(C2A-01A)2 

- > 

(C2A-01A*)2 

0.0121 

(C3A-C4A)2 

- > 

(C3A-C4A*)2 

0.0064 

(01 B-C2B) 2 

- > 

(01 B-C2B*) 2 

0.0093 

(C4A-H6A) 2 

- > 

(C4A-H6A*)2 

0.0066 

(01B-C4A)2 

- > 

(C4A-01B*)2 

0.0036 
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III.  170  NMR  Spectra  of  Oxetanes 

The  170  NMR  spectra  of  several  oxetanes  were  examined  experimentally  by 
Walter  S.  Koski  to  see  if  they  could  be  used  to  calibrate  our  computed 
electrostatic  molecular  potential  contour  (EMPC)  maps  for  the  basicity  of 

oxetane  and  its  derivatives.  The  results  for  oxetane  and  BAMO  show  that  170 
NMR  chemical  shifts  cannot  be  used  directly  to  make  such  calibrations.  The 
reasons  for  this  is  that  two  opposing  effects  contribute  to  the  resultant 
chemical  shift.  One  is  a  diamagnetic  effect,  related  to  the  electron 
density  on  the  9  nucleus,  and  a  second  is  a  paramagnetic  effect,  due  to  the 
presence  of  low  lying  excited  states.  It  is  difficult  to  unscramble  these 
two  contributions  and  get  a  direct  measure  of  the  electron  population  on  the 
oxygen  atom.  However  as  we  point  out  below,  a  detailed  study  of  these  NMR 
shifts  coupled  with  other  experimental  chemicophysical/physiochemical 
measurements  plus  theoretical  interpretation  would  make  valuable 
contributions  to  our  knowledge  of  oxetane  and  its  energetic  derivatives. 

Certain  properties  such  as  basicity  and  electron  density  around  the 
oxetane  play  an  important  role  in  cationic  polymerization  reactions 
involving  oxetane  and  its  energetic  derivatives.  Ab-initio  quantum  chemical 
calculations  and  electrostatic  molecular  potential  contour  (EMPC)  maps 
enable  one  to  get  a  rank  measure  of  the  propensity  of  oxetane  and  its 
derivatives  to  undergo  initiation  of  polymerization.  It  would  be  highly 
desirable  to  determine  experimental 1y  a  property  such  as  electron  density  on 
the  oxygen  atoms  of  some  of  these  compounds.  Such  measurements  would  be  of 
great  value  as  calibration  points  for  the  quantum  chemically  computed 
electrostatic  molecular  potential  contour  maps  to  put  them  on  a  more 
quantitative  absolute  scale.  To  this  end  we  have  made  a  preliminary 

examination  of  using  ^70  NMR  chemical  shifts  to  see  if  they  reflect  the 
electron  density  around  the  oxygen  in  oxetane  and  3,3- 
bis (azidomothyl loxetane  (BAMO).  Neat  samples  of  liquids  were  carefully 
dried  and  sealed  in  10  mm.  tubes  using  D^O  as  an  external  standard.  The 

measurements  were  made  on  a  Varian  400  XL  spectrometer. 

Intuitively  it  would  be  expected  on  the  basis  of  inductive  effects  that 
the  oxygen  in  oxetane  would  have  a  higher  electron  density  and  thus  would  be 
more  shielded  than  the  oxygen  in  BAMO.  This  is  confirmed  by  the  calculated 
electrostatic  molecular  potential  contour  (EMPC)  maps  and  the  calculated 

volumes  enclosed  within  corresponding  isopotential  EMPC  contours.  *70  NMR 
chemical  shifts  in  oxetane  and  BAMO  indicate  that  the  opposite  is  true  for 

the  shielding.  ^70  in  oxetane  falls  10  ppm  lower  than  *70  in  Do0  whereas 

c. 

1 7 

;0  in  BAMO  lies  31.4  ppm  on  the  high  field  side  of  D^O,  indicating  that  the 

0  in  BAMO  is  more  shielded  than  the  0  in  unsubstituted  oxetane.  The  reason 
for  this  is  that  there  are  two  competing  effects  contributing  to  the 
position  of  the  nuclear  magnetic  resonance.  The  first  is  a  diamagnetic 
shielding  effect  caused  by  the  motion  of  electrons  around  the  nucleus  in 
question  (related  to  the  electron  density  on  the  0  nucleus)  and  the  second 
is  a  paramagnetic  term  (due  to  the  presence  of  low  lying  excited  states). 

The  sign  of  the  paramagnetic  term  is  opposite  to  that  of  the  diamagnetic 


effect.  The  effects  of  these  two  terms  are  difficult  to  disentangle 
quantitatively.  This  phenomenon  is  exemplified  by  the  studies  of  Beraldin 

et  al.^  on  the  relationships  between  NMR  shifts  and  atomic  charge  on  ^0  in 
ethers.  Their  results  show  that  in  a  series  of  aliphatic  ethers,  any  gain 

in  electronic  charge  at  the  oxygen  atom  is  accompanied  by  a  downfield 
NMR  shift.  This  is  at  variance  with  the  frequently  held  view  that  any 
increase  in  local  electron  population  should  increase  shielding  of  the 
nucleus  in  question  and  hence  produce  an  upfield  shift. 

2  17 

Similarly  Iwamura  et  al  .  studied  the  NMR  0  shifts  in  twenty-one 
oxiranes  and  found  it  necessary  to  bring  in  diamagnetic,  paramagnetic  and 
other  effects  to  rationalize  the  observed  results. 

However,  a  detailed  and  NMR  study  of  substituted  oxetanes 
coupled  with  other  chemicophysical/physicochemical  measurements  would  give 
us  a  great  deal  of  insight  into  the  electronic  properties  of  the  series  of 
energetic  oxetane  derivatives  of  interest  to  this  ONR  program. 


1.  M.  T.  Beraldin,  E.  Vauthier  and  S.  Flizar.  Can.  J.  Chem.  60,  106  (1982) 

2.  H.  Iwamura,  T.  Sugawara,  Y.  Klawada,  K.  Tori,  R.  Muneyuki  and  R.  Noyori 
Tetrahedron  Letters,  36,  3449  (1979). 
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IV.  Ab-initio  Atom-Class  -  Atom-Class  Potential  Functions  and  Independent 
Confirmations  of  the  Validity  of  Our  Ab-Initio  MODPOT/VRDDO  Energy 
Partitioned  Method  For  Intermolecular  Interactions. 

Over  the  past  several  years  we  had  developed  and  implemented  methods 
for  deriving  ab-initio  atom-class  -  atom-class  potential  functions  from  ab- 

initio  energy  partitioned  SCF  calculations  plus  dispersion  energy.^  The  SCF 
calculations  were  done  with  our  INTER-MOLASYS  program.  That  program  is 
completely  general  for  Gaussian  integrals  either  all-electron  or  with  the 
MODPOT  option  of  using  ab-initio  effective  core  model  potentials.  Our 
earlier  research  on  intermolecular  interactions  calculated  with  our  ab- 
initio  MODPOT  well  balanced  minimal  basis  set  showed  these  basis  sets  have 
comparatively  little  BSSE  (basis  set  superposition  error)  compared  to  all¬ 
electron  calculations.  The  reason  for  this  is  that  while  even  fairly  large 
all-electron  basis  sets  may  describe  the  inner  shells  reasonably  well,  the 
small  deficiency  in  the  inner  shell  wave  function  is  very  heavily  energy 
weighted  and  can  contribute  significant  energy  to  the  BSSE.  On  the  other 
hand,  our  ab-initio  MODPOT  minimal  all-electron  basis  set  had  been  optimized 
to  reproduce  molecular  charge  distributions  resulting  from  larger  well 
balanced  basis  sets.  Then  our  ab-initio  MODPOT  parameter  was  optimized  to 
reproduce  the  valence  electron  molecular  orbital  energies  and  charge 
distributions  of  that  all-electron  well  balanced  minimal  basis  set.  The 
optimization  of  the  MODPCT  basis  set  for  reproducing  molecular  charge 
distributions  of  larger  basis  sets  makes  that  basis  set  more  reliable  than 
other  customarily  used  minimal  atomic  basis  sets  for  calculating 
intermolecular  interaction  energies  of  compounds,  such  as  nitrocompounds  and 
other  compounds  with  electronegative  groups  (such  as  azido.  nitrate,  and 
nitramino)  with  large  intramolecular  charge  redistribution,  since  the 

multipole  part  of  the  electrostatic  component  E^^yp  of  intermolecular 

interaction  energies  is  customarily  the  most  significant  term. 

As  a  minor  effort  of  this  past  year's  ONR  research  (but  with 
significant  favorable  implications  for  our  ONR  research),  we  recently  had 
the  opportunity  to  validate  the  reliability  both  of  the  SCF  intermolecular 
interaction  calculations  using  the  ab-initio  MODPOT  basis  set  and  of  the 
semi-theoretical  expressions  we  had  been  using  for  dispersion  energy  by 
comparing  our  results  on  the  nitromethane  dimer  with  those  of  recent  larger 

2 

basis  set  [4,3, 1/2,1]  all-electron  calculations  on  the  nitromethane  dimer  , 

[SCF,  SCF/SDQ-MBPT(4) ,  SCF/dispersion  energy  Our  ab-initio  MODPOT 

SCF  minimum  energy  occured  at  the  same  intermolecular  distance  as  the  SCF 
minimum  of  the  all-electron  larger  basis  set.  Our  ab-initio  MODPOT  SCF 
interaction  energy  corrected  for  BSSE  (-3.10  kcaVmol  at  2.5  A)  was  only  0.4 
kcal  above  the  SCF  interaction  energy  of  the  larger  basis  set  (-3.57 
keal/mol  corrected  for  BSSE  at  the  minimum  2.5  A). (Figure  7) 
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Furthermore ,  we  also  carried  out  the  ab  initio  MODPOT  <'jnergy- 
partitioned  SCF  calculations  on  the  nitromethane  dimer  to  ascertain  the 
magnitudes  of  the  various  contributions  to  the  interaction  energy.  We  then 

calculated  the  value  of  the  C^MTP  l,sin0  the  multipoles  calculated  for 

f  1 1 SfF 

nitromethane  at  the  SCF  level  [EyL  Myp  with  the  all-electron  larger  basis 

set  used  in  reference  2j.  This  was  done  to  enable  us  to  see  the  genesis  of 
the  small  differences  in  F.<-yp  at  the  SCF  minimum  between  our  ab-initio 

MODPOT  SCF  calculations  and  the  larger  all-electron  SCF  calculations  .  Our 
previous  experience  had  indicated  that  the  major  difference  in 
intermol ecul ar  interaction  energies  between  our  ab-initio  minimal  MOORCT  SCF 
calculations  and  larger  basis  set  calculations  lay  in  a  small  difference  in 

Eel^MTP  ca1cu1d^ed  with  the  different  basis  sets.  Our  present  results  on 

the  nitromethane  dimer  confirmed  this  hyoothesis.  At  the  SCF  minimum 
,  f 1) SCF 

(2.50A)  the  MTP  we  calculated  for  the  nitromethane  dimer  with  the  ab- 

initio  MODPOT  basis  set  was  -4.58  kcal/mole;  the  we  calculated  with 

the  larger  all-electron  [4, 3, 1/2,1]  basis  set  was  -5.09  kcal/mole.  This 

difference,  at  the  SCF  minimum  (2.50  .',)  in  (0.51  kcal/mole)  is  very 

close  to  the  difference  in  SCF  interaction  energies  (0.44  kcal/mole)  at  the 
SCF  minimum.  This  correction  is  certainly  expected  to  be  valid  for  R  >  R  . 

For  the  i ntermol ecul ar  separation  of  2.5  A,  2.7  A,  and  3.0  A  the  correction 

for  the  difference  in  E^MTP  (calcul3teci  explicitly  at  each  point)  brings 

our  total  calculated  interaction  energies  within  0.1  kcal  of  the  AE^y 

values  for  the  larger  basis  set.  Our  ab-initio  MODPOT  SCF  method  (which, 
depending  on  the  nr'Vcules  involved,  is  an  order  to  orders  of  magnitude  more 
rapid  than  al 1-ele- :ron  larger  basis  set  SCF  calculations)  could,  therefore 
be  used  to  calculate  intermoleculur  interaction  energies  to  locate  the  SCF 
minimum.  Our  results  further  suggest  that  it  appears  feasible  to  correct 
for  the  small  differences  in  the  absolute  intermolecular  interaction 
energies  between  the  ab-initio  MODPOT  SCF  results  and  those  of  all-electron 

larger  basis  set  SCF  calculations  by  the  difference  in  the  E^^yp 

calculations  with  the  different  basis  sets. 

In  addition,  the  semi-theoretical  method  we  have  been  using  to  estimate 

dispersion  energies  proved  to  give  values  very  close  to  (variation- 

perturbation )  calculated  by  Cole  et  al .  using  a  larger  all-electron  basis 
set.  (Figure  R) 


Hgure  8 


Nitromethane  Dimer  Dispersion  Energy  [kcal/mol] 

SEMI-THEORETICAL  ATOM-ATOM  POTENTIALS 
Compared  to  E“p  (variation-perturbation  [4,3,1/2,13) 


[Kcal/mol] 
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We  also  explored  another  aspect.  Cole,  et  al .  had  calculated 
the  ^5qq_^bpt(4)  interaction  energy  of  the  nitromethane  dimer  as  the 

difference  between  esdq_MBPT(4) (nitromethane  dimer)  -  -ESDQ_MBpT(4) 

(nitromethane  monomer).  This  e$oq_mBpt  intermolecular  interaction 

energy  includes  both  intermolecular  correlation  energy  and  the  correlation 
contribution  to  the  intermolecular  interaction  energy  from  correlated 

i ntramolecul ar  (monomer)  wave  functions.  We  calculated  ebl^MTP 

for  the  nitromethane  dimer  using  the  correlated  multipoles  from  nitrome¬ 
thane  monomer  which  we  computed  by  carrying  out  correlated  SDQ-MBPT(4) 
calculations  for  nitromethane  monomer  with  the  larger  all-electron  basis  set 
of  Cole  et  al .  excluding  excitations  from  inner  shells. 

CORR 

^INTRA  E|_(correlation  correction  to  intermolecular  interaction 

electrostatic  term  from  multipoles  due  to  use  of  correlated  intramol ecular 
monomer  wave  functions) 


p ( 1 ) SDQ-MDPT (4) 
“  EL,MTP 

basis  set) 


(Cole ,  et  al .  basis  set) 


c(l)SCF 

tEL,MTP 


(Cole ,  et  al . 


»  -4.087  kcal/mole  -(-5.086  kcal/mole) 


=  0.999  kcal/mole  (at  nitromethane  dimer  equilibrium  geometry) 

CORR 

-  AEjujyp^correlation  contribution  to  intermol ecul ar  interaction  energy 

from  intermol ocular  correlation  contributions  due  to  correlated 
monomer  wave  functions) 


-  ( intermolecular  correlation  effects  plus  the  correlation 

contribution  to  i ntermol ecul ar  interaction  energy  from  correlated 
monomer  wave  functions  (at  the  SDQ-MBPT(4)  level)] 


-  E 


(2) 

DISP 


(intermolecular  correlation  effects  (dispersion)} 


*  -0.76  kcal/mole  (at  nitromethane  equilibrium  dimer)  [Cole,et  al.] 

-  1.65  kcal/mole  (at  nitromethane  equilibrium  dimer)  [Cole,  et 
al.] 


-  0.89  kcal/mole. 


This  correction  should  be  valid  for  P  >  R 

eq 


t .  , .  .  .  ...  r ( 1 ) SDQ-MBPT (4)  ,  c(l)SCF 

These  results  suggest  how  one  might  use  MTP  and  ^EL  MTP 

results  to  separate  i ntermol ecul ar  correlation  contributions  to 
intermol ecular  interaction  energies  and  correlation  contributions  to 
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I 
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intermol ecul ar  interactions  energies  which  arise  from  use  of  correlated 
monomer  wave  functions. 


"Crystal  Structure  Studies  Using  Ab-Initio  Potential  Funct.ons  From 
Partitioned  Ab-Initio  MODPCT/VRDDO  SCF  Energy  Calculations.  I.  and 

CO,,  Test  Cases.  II.  Nitromethane W.  Sokalski,  S.  Roszak,  P.  C. 

Hariharan,  W.  S.  Koski,  Joyce  J.  Kaufman,  A.  H.  Lowrey,  and  R.  S. 

Miller.  Int.  J.  Quantum  Chem.,  S17,  375-391  (1983). 

S.  J.  Cole,  K.  Szalewicz,  G.  0.  Purvis,  III  and  R.  J.  Bartlett,  J.  Chem. 
Phys.  84,  6883  (1986) 
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V.  Geometry  Optimizations  at  the  MC-SCF/CI  Level 

We  explored  a  number  of  the  features  of  the  GAMESS  program  which  has 
the  capability  of  optimizing  geometries  by  a  derivative  method  at  the  MC- 
SCF/CI  and  CASSCF  levels  as  well  as  at  the  SCF  level.  As  prototype  problems 
we  have  successfully  optimized  the  geometries  of  oxetane,  protonated  oxetane 
and  nitromethane  at  the  RHF  SCF  level  and  we  have  successfully  optimized  the 
geometry  of  nitromethane  at  the  MC-SCF  level.  (Nitromethane  was  an 
excellent  test  case  since  Dr.  Kaufman  had  previously  used  the  GAMESS  program 
for  MC-SCF  type  calculations  on  the  >C  -  NO^  decomposition  pathway  of 

nitromethane.)  We  will  be  carrying  out  further  studies  to  investigate  the 
effect  of  basis  set  size,  balance,  polarization  functions,  etc.  We  will 
also  be  exploring  and  using  more  of  the  new  options  in  the  GAMESS  program 
for  our  studies  on  the  energetic  oxetanes. 

It  will,  of  course,  always  still  be  necessary  first  to  sample  through 
space  by  carrying  out  the  quantum  chemical  calculations  explicitly  at 
various  geometries  to  find  the  locations  of  the  various  local  maximia  and 
minima  since  geometry  optimizations  by  derivative  methods  lead  only  to  the 
local  minimum  in  the  region  selected  as  the  starting  geometry.  This  is 
where  our  own  ab-initio  MODPOT/VRDDO/MERGE  methods  are  so  useful. 

We  are  also  initiating  derivative  calculations  for  geometry 
optimizations  of  other  substituted  oxetanes  and  substituted  protonated 
oxetane  we  have  investigated  as  a  prelude  to  MRD-CI  calculations  for  their 
interactions.  For  some  of  the  larger  substituted  oxetanes  we  are  having  to 
enlarge  or  modify  the  techniques  in  GAMESS. 

We  are  also  investigating  carrying  out  ab-initio  derivative 
calculations  using  the  GAMESS  program  for  the  most  probable  conformations  of 
the  product  of  the  propagation  step  in  the  cationic  polymerization  of 
oxetane  by  protonated  oxetane.  This  is  a  more  difficult  problem.  In  accord 
with  our  previous  experience  because  of  multiple  maxima  and  minima  one  must 
start  with  a  geometry  in  the  region  of  each  minimum.  Also  if  one  is  too  far 
from  a  minimum,  the  derivative  method  may  not  converge  or  may  converge  to  a 
saddle  point.  We  use  our  ab-initio  MODPOT/VRDDO/MERGE  SCF  techniques  to 
identify  reasonable  starting  conformations 
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VI.  POLY-CRYST 

Commensurate  with  the  ONR  priorities  expressed  to  us  by  our  ONR 
Contract  Monitor  to  devote  our  major  attention  to  the  problem  of  the 
cationic  polymerization  of  oxetanes  ,  we  only  devoted  very  minimal  effort  to 
further  development  and  testing  of  the  POLY-CRYST  program  (for  ab-initio 
calculations  of  crystals  and  polymers). 

For  efficiency  we  initiated  changes  in  the  way  the  integrals  are  picked 
up.  We  also  ran  a  few  test  calculations.  The  stragegy  for  testing  is  to 
calculate  the  energy  of  the  isolated  molecule  and  the  cluster  of  molecules 
in  the  unit  reference  cell  with  our  ab-initio  molecular  (MOLASYS)  or 
intermolecular  programs  ( INTER-MOLASYS) .  Then  we  repeat  the  calculation  on 

the  unit  reference  cell  with  POLY-CRYST  with  i<  =  0.  The  energy  of  the 
cluster  of  molecules  in  the  unit  reference  cell  (by  INTER-MOLASYS)  and  that 
of  the  unit  reference  cell  by  POLY-CRYST  must  agree.  We  then  run  the  POLY- 
CRYST  calculations  for  a  number  of  cells. 

The  POLY-CRYST  program  has  great  potential  for  understanding  properties 
of  energetic  molecular  crystals.  We  hope  that  in  the  future  ONR  priorities 
will  be  such  that  we  can  devote  more  attention  to  this  research  area. 
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VII.  Lectures  Presented  and  Publications  on  This  ONR  Research 

Presentations  given  and  or  scheduled  and  papers  published  and/or 
submitted  during  the  fiscal  year 

A.  Presentations  Given  Dr.  Joyce  J.  Kaufman 

★ 

1.  Already  presented  (  denotes  invited  lecture) 

a.  National  and  International  Meetings 

*  "Ab-Initio  MRD-CI  Calculations  of  the  >C  -  N0o 

<1 

Decomposition  Pathway  of  Nitrobenzene,"  Sanibel 
Internat ional  Symposia  on  Quantum  Chemistry,  Solid  State 
Theory,  Many-Gody  Phenomena  and  Computational  Quantum 
Chemistry,  Marineland,  Florida,  March  1986 

*  "Ab-Initio  Multireference  Double  Excitation  Configuration 
Interaction  Calculations  Based  on  Localized  Orbitals  for 
Molecular  Decompositions  and  Reactions  for  Large  Systems," 
Sanibel  International  Symposia  on  Quantum  Chemistry,  Solid 
State  Theory,  Many-Body  Phenomena  and  Computational  Quantum 
Chemistry,  Marineland,  Florida,  March  1986. 

"More  New  Desirable  Computational  Strategies  for  Ab-Initio 
Calculations  on  Large  Molecules,  Clusters,  Solids  and 
Crystals,"  Sanibel  International  Symposia  Snowbel 
Satellite  on  Interface  Between  Electronic  Structures  and 
Dynamics,  Snowbird,  Utah,  April  1986. 

*  "Ab-Initio  MRD-CI  Calculations  on  the  >C  -  N0£ 

Decomposition  Pathway  of  Nitrobenzene,"  plenary  lecture 
presented  at  the  Symposium  on  Computational  and 
Mathematical  Chemistry,  Can.  Inst.  Chem.  National  Meeting, 
Saskatoon,  Canada,  June  1986. 

*  "Ab-Initio  Potentials  for  Crystals  and  Ab-Initio  Crystal 
Orbitals,"  special  invited  lecture  presented  at  the 
International  Symposium  on  Molecules  in  Physics,  Chemistry 
and  Biology,  dedicated  to  Professor  R.  Deudei  ,  Paris, 
France,  June  1986. 

*  "Ab-Initio  Localized  Orbital  MRD-CI  Calculations  on  the 
Mechanism  of  Cationic  Polymerization  of  Oxetane  By 
Protonated  Oxetane,"  IUPAC  5th  Internat i onal  Symposium  on 
Ring-Opening  Polymerization,  Blois,  France,  June  1986. 

*  "Ab-Initio  Multirefcrence  Double  Excitation  Conf iguration 
Interaction  Calculations  Based  on  Localized  Orbitals  for 
Molecular  Decompositions  and  Reactions  for  Large  Systems,: 
American  Chemical  Society  National  Meeting,  Anaheim, 

Calif.,  September  1986. 
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b.  Other  Research  Institutions 

*  "Ab-Initio  Calculations  on  Large  Molecules,  Crystals  and 
Solids  Using  Desirable  Computational  Strategies,"  Sandia 
Laboratories ,  Albuquerque,  New  Mexico,  May  1986. 

*  "Ab-Initio  Calculations  on  Large  Molecules,  Crystals  and 
Solids  Using  Desirable  Computational  Strategies,"  Los 
Alamos  National  Laboratory,  New  Mexico,  May  1986. 

2.  To  be  Presented 

a.  National  and  International  Meetings 

*  "Ab-Initio  MRD-CI  Calculations  for  the  Prediction  of 
Cationic  Polymerization  of  Oxetanes  Based  on 
Localized  Orbitals,"  Sanibel  International  Symposia  on 
Atomic,  Molecular  and  Solid  State  Theory,  Many-Body 
Phenomena  and  Computational  Methods,  Marineland,  Florida, 
March  1987. 

*  "Comparison  of  Ab-Initio  MODPOT  and  Ab-Initio  Energy 
Partitioned  Potential  Functions  for  Nitromethane  Dimer 
Against  Large  Basis  Set  Calculations,"  Sanibel 
International  Symposia  on  Atomic,  Molecular  and  Solid  State 
Theory,  Many-Body  Phenomena  and  Computational  Methods, 
Marineland,  Florida,  March  1987. 

*  "Ab-Initio  Calculations  on  Large  Molecules  and  Solids  Using 
Desirable  Computational  Strategies,"  Vlllth  International 
Conference  on  Computers  in  Chemical  Research  and  Education, 
Beijing,  China,  June  1937. 

b.  Other  Research  Organi zations 

*  "Ab-Initio  Quantum  Chemical  Calculations  on  Large  Molecular 
Systems  and  Crystals",  Maryland  Section  American  Chemical 
Society,  February  1987. 

c.  At  Department  of  Defense  Meetings 


"Ab-Initio  Multireference  Determinant  Configuration 
Interaction  Calculations  for  the  Propagation  Step  in 
Cationic  Polymerization  of  Oxetanes,"  ONR  Workshop  on 
Polymeric  Energetic  Materials  Synthesis,  Great  Oak, 
Maryland,  October  1986. 


"Ab-Initio  MRD-CI  (Mul ti ref erence  Double  Excitation 
Configuration  Interaction)  and  Ab-Initio  Intermolecul ar 
Interaction  Calculations  on  Energetic  Nitrocompounds,"  ONR 
Workshop  on  Crystalline  Energetic  Materials  Synthesis, 
Great  Oak,  Maryland,  November  1986. 


.'iVaW, 


Papers  Published  on  This  ONR  Research,  October  1985-September  1986 

1.  Already  Published 

a.  “Ab-Initio  Multireference  Determinant  Configuration 
Interaction  (MRD-CI)  and  CASSCF  Calculations  on  Energetic 
Compounds,"  Joyce  J.  Kaufman,  P.  C.  Hariharan,  S.  Roszak, 

C.  Chabalowski,  M.  van  Hemert,  M.  Hotokka  and  R.  J. 

Buenker.  In  CHEMICAL  REACTION  DYNAMICS,  NATO  Advanced  Study 
Institute,  Editors  P.  M.  Rentzepis  and  C.  Capellos,  D. 
Reidel  Publishers,  Dordrecht ,  Netherlands,  1986,  pp.  289- 
310. 

b.  "Crystal  Structures  of  Energetic  Compounds:  Ab-Initio 
Potential  Functions  and  Ab-Initio  Crystal  Orbitals,"  Joyce 
J.  Kaufman,  P.  C.  Hariharan,  S.  Roszak,  J.  M.  Blaisdell  and 
A.  H.  Lowrey,  and  R.  S.  Miller.  In  CHEMICAL  REACTION 
DYNAMICS,  NATO  Advanced  Study  Institute,  Editors,  P.  M. 
Rentzepis  and  C.  Capellos,  D.  Reidel  Publishers, 

Dordrecht,  Netherlands,  1986,  pp.  311-326. 

c.  "POLY-CRYST  -  A  Program  for  Ab-Initio  Crystal  Orbitals  and 
Polymer  Orbitals,"  J.  J.  Blaisdell,  W.  A.  Sokalski,  P.  C. 
Hariharan  and  Joyce  J.  Kaufman.  J.  Non-Crystalline  Solids, 
75,  319-326  (1985). 

d.  "Procedure  Supplementing  SCF  Interaction  Energies  By 
Dispersion  Term  Evaluated  in  Dimer  Basis  Within  Variation 
Perturbation  Approach,"  S.  Roszak,  W.  A.  Sokalski,  P.  C. 
Hariharan  and  Joyce  J.  Kaufman.  Theo.  Chim.  Acta.,  70,  81- 
88  (1986). 

2.  Accepted  For  Publication  and  In  Press 

a.  "Ab-Initio  Localized  Orbital  MRD-CI  Calculations  on  the 
Mechanism  of  Cationic  Polymerization  of  Oxetane  By 
Protonated  Oxetane,"  Joyce  J.  Kaufman,  P.  C.  Hariharan,  S. 
Roszak  and  M.  Hotokka.  An  invited  lecture  presented  at  the 
IUPAC  5th  International  Symposium  on  Ring-Opening 
Polymerization,  Blois,  France,  June  1986.  In  Press, 
Makromolecular  Chemie,  special  Symposium  Issue. 

b.  "Ab-Initio  MRD-CI  Calculations  on  the  >C-NO?Decomposition 

Pathway  of  Nitrobenzene,"  Joyce  J.  Kaufman,  P.  C. 

Hariharan,  S.  Roszak  and  M.  van  Hemert.  An  invited  plenary 
lecture  presented  at  the  Symposium  on  Computational  and 
Mathematical  Chemistry,  Can.  Inst.  Chem.  National  Meeting, 
Saskatoon,  Canada,  June  1986.  In  press,  J.  Comp.  Chem., 
Symposium  Issue. 

c  "Symposium  Note:  More  New  Desirable  Computational 

Strategies  for  Ab-Initio  Calculations  on  Large  Molecules, 
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Clusters,  Solids  and  Crystals,"  Joyce  J.  Kaufman.  In 
press,  Int.  J.  Quantum  Chem.,  Symposium  Issue. 

d.  "Nonempirical  Atom-Atom  Potentials  for  Main  Components  of 
Intermolecular  Interaction  Energy,"  W.  A.  Sokalski,  A.  H. 
Lowrey,  S.  Roszak,  V.  Lewchenko,  J.  M.  Blaisdell,  P.  C. 
Hariharan  and  Joyce  J.  Kaufman.  In  press,  J.  Comp.  Chem. 
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VIII.  Project  Personnel 

Joyce  J.  Kaufman,  Ph.D. 

Principal  Investigator 


P.  C.  Hariharan,  Ph.D. 

Research  Scientist 

Overall  responsibility  for  implementing  new  program  developments 
and  conversion  to  CRAY  computers. 

Quantum  chemical  calculations  on  energetic  polymers,  MRD-CI, 
GAMESS  and  POLY-CRYST  calculations 


Philip  B.  Keegstra,  Ph.D. 

Post  Doctoral 

MRD-CI  calculations  on  cationic  polymerization  of  oxetanes,  test 
calculations  on  POLY-CRYST. 


Robert  J.  3uenker,  Ph.D. 

Visiting  Col laborator  (March-April  1986) 

Initiated  expansion  of  MRD-CI  program  to  include  more  rrference 
configurations,  more  roots  and  a  post  MRD-CI  population  analysis  for 
NRL  CRAY  XMP-12  (COS). 

(Permanent  Address:  University  of  Wuppertal,  West  Germany) 


Marc  van  Hemert,  Ph.D. 

Visiting  Collaborator  (June-July  1986) 

Continued  expansion  of  MRD-CI  program  and  converted  it  to  CRAY 
XMP-48  (CTSS) 

(Permanent  address:  University  of  Leiden,  Netherlands) 


S.  Roszak  Ph.D. 

Visiting  Associate  Research  Scientist  (July  1985  -  January  1986) 

Test  MRD-CI  calculations  to  validate  use  of  localized  orbitals 
MRD-CI  calculations  for  cationic  polymerization  -  opening  the  oxetane 
ring  (neutral  and  protonated). 


W.  A.  Soka 1  ski ,  Ph .D. 


Visiting  Associate  Research  Scientist  (August-September  1986) 

Program  developments  for  correlated  wave  functions.  Test  of  our 
MODPOT  and  additional  methods  against  large  basis  set  calculations  for 
nitromethane  dimer. 
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